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SUMMARY 


Generalized  Youden  Designs:  Construction  and  Tables 

toy 

Arlene  Ash 

University  of  Illinois  at  Chicago 
and 

Boston  University 

Generalized  Youden  Designs  are  generalizations  of  the  class 
of  two-way  balanced  block  designs  which  include  Latin  squares 
and  Youden  squares.  They  are  used  for  the  same  purposes  and  in 
the  same  way  that  these  classical  designs  are  used,  and  satisfy 
most  of  the  most  common  criteria  of  design  optimality. 

We  explicitly  display  or  give  detailed  instructions  for 
constructing  all  these  designs  within  a practical  range:  when 
v,  the  number  of  treatments,  is  £ 25;  and  bj^  and  b2,  the 
dimensions  of  the  design  array,  are  each  ^ 50. 

* ••  • 
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2 

GYD(v,b^,b2)  and  GYD(v' jb^.bl,)  are  in  the  same  family  if  v = v'  and  bA 
is  congruent  to  b^  mod  v,  for  i = 1,2. 

3 

The  rows  of  a Bq  ■ GYD(22,22,7)  would  be  a BIB(22,22,7)  which  cannot  be 
constructed.  See  for  example  (Collens,  1976).  A design  with  twice  as  many 
blocks,  the  smallest  possible  in  this  family,  is  listed  instead. 

4 

No  small  design  In  this  family  exists,  because  of  the  non-cxis ter.ee  of 
required  BIB(15 ,21 ,5) . See  for  example  (Collens,  1976). 


by 

Arlene  Ash 


University  of  Illinois  at  Chicago 
and 

Boston  University 

1.  Introduction . A Generalized  Youden  Design  (GYD)  is  a * b2  array 
on  the  treatments  denoted  GYD(v,b1 ,b2) , whose  induced  row  and 

column  designs  satisfy  certain  technical  conditions  of  balance.  It  is  a 
natural  generalization  of  several  well-known  highly  balanced  two-way  design 
structures,  such  as  the  Latin  square  and  the  Youden  square.  Indeed,  when 
bj  ■ v and  b2  S v,  a GYD(v,b^,b2)  is  exactly  one  of  these  designs.  The  GYD 
definition  extends  the  collection  of  design  settings  in  which  two-way  bal- 
ance may  be  achieved,  to  include  values  of  b^  and  b2  which  may  be  greater 
than  v. 


GYD's  are  uccd  for  comparing  the  performance  of  v treatments  in  the 
presence  of  two-way  heterogeneity  (in,  respectively,  b1  and  b2  blocks)  under 
the  standard  linear  model  with  no  interactions.  For  such  purposes  GYD's  are 
optimal  under  several  of  the  most  common  criteria  of  optimality  [6].  This 
paper  makes  these  statistically  good  designs  available  by  listing  (or  giving 
explicit  directions  for  constructing)  those  GYD's  within  a practical  range: 
v S 25;  b1,b2  s 50. 

Associated  with  each  design  is  a v x v matrix  called  the  information 
matrix  of  the  design  (for  the  treatment  means).  Usually  the  goodness  of 
a design  is  measured  as  some  function  of  this  matrix.  Since  non-isomorphic 
GYD's  with  the  same  parameters  (v,bj,b2)  have  identical  information  matri- 
ces, we  indicate  only  one  possible  construction  for  each  parameter  set  where 
construction  can  be  achieved.  In  the  range  specified  there  is  only  one 


I 


1.2 


parameter  set  (v  = 25,  = \>2  = AO)  where  neither  nonexistence  of  the  GYD, 

nor  any  construction  is  known.  We  have,  however,  provided  a 40  x AO  design, 
P18,  the  combined  set  of  whose  rows  and  columns  form  a balanced  80-block 
design.  This  makes  it  a Pseudo  Youden  Design  (PYD)  as  defined  by  Cheng 
(1978)  and  gives  it  the  same  information  matrix  and  optimality  properties 
as  the  missing  GYD. 

Previous  work  by  Seiden,  Ruiz  and  Wu  (197A,  1978)  has  provided  explicit 
information  for  constructing  many  designs  where  v is  a power  of  a prime. 
Theorems  by  Kiefer  (1975)  combined  with  lists  of  BIB's  (Hanani,  1975  or 
Collens,  1976  for  example)  guarantee  the  existence  of  some  other  designs 
in  our  listing.  Designs  P18  and  N3,  6,  10,  11,  12,  1A  and  16  are  apparently 
new. 

Special  acknowledgements  are  due  to  Ching-Shiu  Cheng  who  pointed  out 
that  a BIB (25, 80, 15)  could  be  substituted  for  two  BIB (25  ,40 , 15) 1 s in  making 
a PYD  surrogate  ror  GYD (25 , AO , AO) , and  to  Haim  Hanani  who  provided  the  eighty 
block  design  which  was  used  in  making  that  construction.  Thanks  also  to 
R.  J.  Collens  who  helped  locate  a needed  BIB  design  from  the  updated  listing 
he  maintains  on  the  University  of  Manitoba  computer. 

It  still  seems  of  some  theoretical  interest  to  know  whether  or  not 
any  GYD(25 ,40 ,40)  exists. 


2. 


Definitions  and  notation.  Throughout  this  paper  we  will  let 


l>i  = m^v  + cif  0 <;  Ci  < v,  1=  1,2. 

Definition.  A GYD(v,b^,b2)  is  a b|  k b2  array  with  entries  from 
V = {l,...,v}  such  that 

(a)  Every  element  of  V occurs  either  m^  or  m^  + 1 times  in  each  column 
and  either  m2  or  ®2  + 1 times  in  each  row;  and 

(b)  Every  two  distinct  elements  of  V occur  together  in  the  same  row  the 
same  number  of  times;  and  similarly  for  columns. 

That  is,  the  rows  of  a GYD  form  a so-called  balanced  block  design, 

BBD(v,b=bj ,k=b2) ; and  the  columns  form  a BBD (v ,b=b2 ,k=b ^) . The  existence 
of  BBD's  with  the  given  parameters  is  equivalent  to  the  existence  of  bal- 
anced incomplete  block  designs  BIB(v,b=b^ ,k=c2)  and  BIB(v,b=b2 ,k=c^) . Well- 
known  parameter  conditions  are  necessary  for  the  existence  of  BIB's,  and 
consequently,  for  a fixed  (v,Cj,c2),  GYDCv.b^m^v+c^  ,b2=m2v+c2)  can  only 
exist  when  v divides  c^c2  and  the  b^,  in  addition  to  being  congruent  to 
c^  mod  v,  satisfy  these  BIB-related  conditions. 


Construction.  GYD  construction  is  trivial  when  v divides  both  bj  and 


b2i and  essentially  equivalent  to  the  construction  of  a related  BIB  when  v 
divides  just  one  of  or  b2  but  not  both.  It  is  new  and  interesting  when 
neither  b^  nor  b2  is  a multiple  of  v.  We  distinguish  these  three  design 
settings  as,  respectively,  the  doubly  regular,  (simply)  regular  and  non- 
regular cases. 


3.1  The  doubly  regular  case.  GYD (v,b^=m^v,b 2=11^)  can  be  constructed  by 
juxtaposing  any  m^  by  m2  collection  of  Latin  squares.  For  example,  let 
the  (i,j)th  entry  of  the  array  be  v if  i + j = 1 mod  v and  (i  + j - l)mod  v 
otherwise. 


3.2  The  (simply)  regular  case.  Transposing  the  design  if  necessary  we 
may  assume  that  c^  = 0.  We  may  construct  GYD(v  ,b^=mjV  ^2=11^+02)  by  ad- 
joining a doubly  regular  GYD^.mjV.ii^v)  to  a regular  GYD^.m^v  ,c^)  . 

For  every  v and  0 < < v,  there  is  a b^  = m^v,  such  that  if  regular 

Bj  = GYD(v,bj=nijV,C2)  exists  then  b^  is  a multiple  of  b^.  If  B^  = GYD(v,bQ,C2) 
exists,  then  can  be  constructed  by  stacking  bj/bp  copies  of  B^. 

= Bq(v,C2)  is  called  a core  design  for  the  (v , C2) -family  of  regular 
designs.  We  list  instructions  for  finding  B^: 


3.2.1.  If  C2  = 1,  Bq  may  be  taken  as  the  (v  * 1)  column  array  C 1 ,2, . . . ,v3 ' , 
where  the  "prime"  indicates  transpose. 


3.2.2.  If  C2  a 2,  Bq  will  be  a v(v  - 1)  x 2 array  for  v even,  and  a 
[v(v  - l)/2]  x 2 array  for  odd  v.  In  either  case,  construct  the  first 
column  of  BQ  by  stacking  copies  of  column  vector  D = [ 1 ,2 , . . . ,v] ' ; and  the 
second  column  by  cycling  the  entries  of  D : i.e. , first  use  [2,3, . . . ,v, 1]  ' , 
then  [3,4 , . . . ,v, 1 , 2] ' , etc.  until  the  second  column  is  filled.  E.g.  for 
v = 5 we  get 


Bq  = Bq(5,2)  = GYD (5 ,10,2)  = 


3.2.3.  If  3 S C2  ^ v/2,  the  construction  of  Bq  is  more  complicated.  For 

v s 25  and  Bq  no  bigger  than  50,  the  design  Eq  will  be  found  in  the  regular 
design  tables  through  the  parameters  (v,C2). 

3.2.4.  For  C2  > v/2,  there  is  a three-step  procedure.  First,  locate  Bq, 
the  core  design  for  the  (v,C2=v-C2)  family,  as  above.  Then,  construct 
coiiip(Bq),  a row  design  complementary  to  Bq  viewed  only  as  a row  design. 

That  is,  put  i e {l,...,v}  in  the  j*"*1  row  of  the  complementary  design  when 

i is  not  in  the  row  of  the  original.  This  new  design  is  a BIB (v ,b ^ , C2)  , 
but  the  array  may  require  rearrangement  of  the  entries  within  rows  so  that 
the  columns  are  also  balanced.  This  can  always  be  done  (Agrawal , 1966), 
and  for  small  designs  it  can  be  done  rapidly  by  eye.  Sometimes,  if  Bq  has 
been  constructed  with  a discernible  column  pattern  this  facilitates  writing  B 


3.3 


for  GYD(7 , 21 , 18)  we  need  GYD(7 ,21 ,4) . Since 


directly  as  a GYD.  E.g 


A 21  x 14  doubly  regular  design,  patched  to  a stack  of  three  copies 


of  the  second  version  of  comp(B') , makes  the  desired  GYD 


If  "eyeball"  methods  fail  to  produce  a rearrangement  of  the  BIB  which 


is  column-balanced  the  appendix  contains  a definitive  algorithm  (called 


SHUFFLE)  for  turning  a BIB  with  b a multiple  of  v into  a GYD 


3.3.  The  Nonregular  Case 


when  v | 0^2 


For  such  a triple  (v,c.,c9)  there  are  integers  B , 6 


such  that  any  GYD(v,b.,b9)  with  b 


sible  dimensions  in  the  (v,c. ,c9)-family ; and  the  b 


for  which  BIB's  in  b 


blocks  of  size  c.  can  exist.  G 


it  exists,  is  called  a core  design  for  the  family.  If,  further,  regular 


family  can  be  constructed.  One  method  is  indicated  in  the  diagram  below 


where  the  hatched  regions  are  doubly  regular  designs 
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Since  nonregular  designs  cannot  readily  be  assembled  from  the  asso- 
ciated row  and  column  designs,  explicit  constructions  for  the  small  core 
designs  are  given  in  4.3.  The  parameters  of  the  associated  designs  are 
indicated  next  to  those  of  in  the  listing  which  precedes  the  construc- 

tions. Most  of  these,  however,  have  too  many  blocks  to  be  included  in  the 


tables  of  regular  designs. 


A. 


Guide  to  the  tables.  The  tables  contain  the  building  blocks  for 


, 


A . 1 


GYD(v,b^,b2) 's  with  v S 25  and  b^.b^  both  S 50.  We  refer 
as  "small".  In  order  to  determine  whether  or  not  a small 
to  find  its  construction  if  it  does,  perform  the  numbered 
until  they  lead  to  an  exit. 


to  such  designs 
GYD  exists  and 
steps  sequentially 


The  exits  are  coded  by 

2 

QD  GYD(v,bpb2)  does  not  exist,  nor  does  any  GYD  in  the  same  family  exist. 

0 GYDCv.bpb^)  with  b^  and  b^  both  < 50  does  not  exist  although  larger 

GYD's  in  the  same  family  will. 

1=S>  Proceed  elsewhere,  as  directed. 

The  procedure  is  definitive  for  v £ 25  and  b^jb^  both  < 50,  although 
a Pseudo  Youden  Design  (see  the  Introduction)  has  been  offered  instead  of 
GYD(25, AO, AO)  in  the  one  place  where  we  have  not  been  able  to  achieve  the 
more  rigid  construction.  Even  if  b^  or  b2  is  bigger  than  50,  if  v < 25  a 
construction  may  result  if  the  GYD  sought  is  in  the  same  family  as  another 
design  which  is  in  the  range  of  the  tables.  For  v > 25,  the  procedure  is 
only  useful  for  determining  when  no  design  can  exist. 


A.l.  Procedure  for  resolving  a GYD(v,b^,b2)  construction  problem. 

1)  Compute  c^,  the  residue  of  b^  mod  v,  for  i = 1,2. 

2)  If  v'fc^C2»  then  no  design  exists,  dl 

3)  If  c^  = Cj  = 0,  the  design  construction  is  doubly  regular  and  can 
be  achieved.  See  3.1. 

2 

GYD(v,bj,b2)  and  GYD(v’ ,bj  ,bp  are  in  the  same  family  if  v = v1  and  b^ 
is  congruent  to  b^  mod  v,  for  i = 1,2. 


1 


w 


4.2 


4)  If  necessary  relabel  and  b ^ so  that  Cj  s c^. 

^1C2  blc2^c2-^ 

5)  If  — — and  are  not  both  integers  then  no  design  exists.  CD 

6)  If  Cj  = 0,  then  the  construction  problem  is  regular.  In  this  case 
do  the  following: 

a)  If  C2  = 1 or  2 the  design  construction  is  trivial.  See  3.2.1 
or  3.2.2,  respectively.  I — \ 

b)  Compute  h2  = g.c.d. (v-l,c2(c2-l)) • 

c)  If  v(—r — -)  > 50,  no  small  design  exists.  □ 

h2 

d)  If  3 s Cj  ^ v/2,  look  for  the  design  in  4.2. 

e)  If  c 2 > v/2,  see  3.2.4. 

b2ci  b2Cl^Cl  “ ^ 

7)  If  — — and  — — j are  not  both  integers,  then  no  design 

exists.  CD 

8)  Compute  g^,  = g.c.d. (v,c^)  and  h^  = g. c. d . (v-1  .c^ (c^-1) ) , for  i = 1,2. 

9)  If  either  V^— ■ — — or  -V^V  — — is  > 50,  no  small  design  exists.  □ 

g2h2  8lhl 

10)  If  v = 25  and  b^  = b2  = 40,  it  is  not  known  whether  a GYD  exists. 

An  equivalent  design  PYD(25,40,40)  is  listed  in  4.3  as  design  P18. 

11)  Look  for  the  design  in  4.3. 


4.2.  Listings  for  regular  designs.  (v  £ 25;  c^  =0;  C2  *0). 

For  GYD(v,bj=m^v,b2=m2V+C2)  to  exist  bj  must  be  a multiple  of 


b0  “ b0(v*c2) 


v(v  - 1) 

.c.d.  (v-l,c2(c2-l))  ’ 


where  g.c.d.  denotes  greatest  common  divisor. 


For  v i 25,  with  3 i Cj  ^ v/2,  if  £ 50,  this  table  displays 
Bq  = GYD(v,bQ, c2) , from  which  G can  be  made  by  stacking  bj/b^  copies  of 


Bq  and 

adjoining 

an  m^v  x n^v  doubly 

regular 

design  (see  3.1).  For  one 

parameter  set  in 

design  with  2bp  = 

this  range,  v = 22 

s 44  blocks  is  given 

and  c2  = 

instead 

7,  Bq  does  not  exist,  and  , 

For 

v s 25  and  c^  > v/2,  Bq  can 

be  constructed  as 

the  complement  of 

a (v,c^ 

) -design, 

where  c^  = v - c2» 

(See  3. 

2.3.) 

For 

any  v and 

1 c2  = 1 or  2,  see  sections 

3.2.1  or 

3.2.2  respectively 

Listing 

Family 

Dimensions  of  B 

Listing 

Family 

Dimensions  of  B 
GYD(v,b0,c2)  U 

Name 

(v,c2) 

GYD(v,b0,c2)  U 

Name 

(v,c2) 

R1 

(6,3) 

(6,30,3) 

R13 

(16,5) 

(16,48,5) 

R2 

(7,3) 

(7,7,3) 

R14 

(16,6) 

(16,16,6) 

R3 

(9,3) 

(9,36,3) 

R15 

(17,8) 

(17,34,8) 

R4 

(9,4) 

(9,18,4) 

R16 

(19,9) 

(19,19,9) 

R5 

(10,3) 

(10,30,3) 

R17 

(21,5) 

(21,21,5) 

R6 

(10,4) 

(10,30,4) 

R18 

(21,6) 

(21,42,6) 

R7 

(11.5) 

(11.11,5) 

R19 

(21,10) 

(21,42,6) 

R8 

(13,3) 

(13,26,3) 

R20 

(22,7) 

(22,44,7)3 

R9 

(13,4) 

(13,13,4) 

R21 

(23,11) 

(23,23,11) 

RIO 

(13,5) 

(13,39,5) 

R22 

(25,4) 

(25,50,4) 

Rll 

(13,6) 

(13,26,6) 

R23 

(25,9) 

(25,25,9) 

R12 

(15,7) 

(15,15,7) 

R24 

(25,12) 

(25,50,12) 

Table  1:  Small,  core  regular  designs. 


The  rows  of  BQ  = GYD(22,22,7)  would  be  a BIB(22,22,7) , which  cannot  be 
constructed.  See  for  example  (Collens,  1976).  A design  with  twice  as 
many  blocks,  the  smallest  possible  in  this  family,  is  listed  instead. 
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4.3.  Listings  for  nonregular  designs  (v  5 25;  0 * c^  < c^). 

For  GYD(v,  “ m^v  + c^,  ° n^v  + Cj)  to  exist  v must  divide  CjC2> 

and  the  b^'s,  in  addition  to  satisfying  congruency  to  c^  mod  v,  must  have, 
for  i = 1,2  and  j = 2,1  respectively,  b^  is  a multiple  of 

b = y . (.y  - -1-) 

i g.c.d.(v.Cj)  g.c .d . (v-1 ,Cj (Cj-1)  ’ 

where  g.c.d.  denotes  greatest  common  divisor. 

Core  designs  Gq  * GYD(v, 8^ , B2)  in  this  family  are  listed  below,  whenever 
these  smallest  possible  constructions  have  v < 25,  8j  and  (32  - 50.  Larger 
designs  in  the  same  family  can  be  constructed  by  making  a patchwork  of  regu- 
lar designs,  using  B°  = GYD(v,bJ ,c2)  and  B°  = GYD(v,b°,Cj)  in  the  construction. 
CSee  3.3.) 

Table  2: 

Small,  core  nonregular  designs. 


Listing 

Name 

(v,clfc2)- 
f amily 

Core  Design,  G 
(v.Bj.Sj) 

Associated 
= (v,b°,c2) 

Designs 

B°  = (v,b°,Cj) 

N1 

(4,2,2) 

(4,6,6) 

(4,12,2) 

(4,12,2) 

N2 

(6,2,3) 

(6,20,15) 

(6,30,3) 

(6,30,2) 

N3 

(6,3,4) 

(6,15,10) 

(6,30,4) 

(6,30,3) 

N4 

(8,2,4) 

(8,42,28) 

(8,56,4) 

(8,56,2) 

N5 

(8,4,4) 

(8,28,28) 

(8,56,4) 

(8,56,4) 

N6 

(8,4,6) 

(8,28,14) 

(8,56,6) 

(8,56,4) 

N7 

(9,3,3) 

(9,12,12) 

(9,36,3) 

(9,36,3) 

N8 

(9,3,6) 

(9,12,24) 

(9,36,6) 

(9,36,3) 

N9 

(9,6,6) 

(9,24,24)  , 

(9,36,6) 

(9,36,6) 

N10 

(10,5,6) 

(10,15,36) 

(10,30,6) 

(10,90,5) 

Nil 

(10,5,8) 

(10,45,18) 

(10,90,8) 

(10,90,5) 

N12 

(12,8,9) 

(12,44,33) 

(12,132,9) 

(12,132,8) 

_ 4 

(15,5,6) 

(15,35,21) 

(15,105,6) 

(15,105,5) 

N14 

(15,5,12) 

(15,35,42) 

(15,105,12) 

(15,105,5) 

N15 

(16,4,4) 

(16,20,20) 

(16,80,4) 

(16,80,4) 

N16 

(21,9,14) 

(21,30,35) 

(21,210,14) 

(21,150,9) 

N17 

(25,5,5) 

(25,30,30) 

(25,150,5) 

(25,150,5) 

P185 

(25,40,40) 

(25,40,40) 

(25,300,15) 

(25,300,15) 

No  small  design  in  this  family  exists,  because  of  the  non-existence  of 
required  BIB(15, 21,5) . See  for  example  (Collens  , 1976). 


The  construction  listed  is  a Pseudo  Youden  Design.  See  the  Introduction. 
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80  PRINT*  FILE  NAMES?  * 

90  INPUT  A $ * D t 
100  FILE  #1JA* 

110  FILE  821B$ 

115  MARGIN  *210 

120  DIM  A (50 *25) *N(50*25> *G<50> 
130  INPUT  4M1V*B*K 
140  MAT  A=ZER(B*K> 

150  MAT  N-ZER ( B * K ) 

160  MAT  INPUT  tllA(B*K) 

170  LET  M=INT<B/V> 


180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

200 

290 

300 

310 

320 

330 

332 

334 

336 

338 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 


FOR 


P=1 

FOR 


TO  U 

J-l 

FOR 


TO  K 

1 = 1 TO  B 

IF  A(I*J)--=P 

THEN  LET  A < I * J) =P+C*V 

N ( P+C*V * J ) --N ( P4-C*V * J)  4 1 
C =MOn<C+l *M) 


LET 

LET 

CONTINUE 

I 


NEXT 
NEXT  J 
IF  OO 

THEN  PRINT  "IMPROPER 
LET  0=1 
CONTINUE 
NEXT  P 

LET  M=1 
PERFORM  CHECK 
DO  WHILE  M=1 

PERFORM  CHECK 
FOR  C=1 
FOR 


* OF  TREATMENTS* *P 


TO  K 

R=1  TO  B 
IF  N ( R * C ) =0 
THEN  LET  Cl=l 

DO  WHILE  N( R * Cl ><=1 
LET  C1=C1 41  . 

IF  C1>K.  THEN  PRINT  ’HELP* 
LOOP 

MAT  G=ZER ( B > 

LET  X=R 
PERFORM  FINDY 

DO  WHILE  N ( Y * C ) = 1 AND  N(Y*C1)>0 
PERFORM  EXCHANGE 
LET  X=Y 
PERFORM  FINDY 
LOOP 

PERFORM  EXCHANGE 
CONTINUE 


520  NEXT  R 

530  NEXT  C 
535  LOOP 


040 

IF  0=0 

545 

PERFORM  UNSCRAMBLE 

j 

550 

THEN 

570 

SCRATCH  12 

575 

PRINT  #2:V» * r 1 JBf " r 1 »K 

580 

FOR  1=1  TO  B 

590 

FOR  J=1  TO  K 

600 

PRINT  I2JA(I r J) i 1 f 

■ ; 

610 

NEXT  J 

620 

PRINT  *2J 

» 

630 

NEXT  I 

640 

ELSE  PRINT  ‘NO  WRITING  DONE  IN 

FILE*  f Bt 

650 

CONTINUE 

670 

680 

DEFINE  FINDY 

690 

LET  1=1 

700 

DO  WHILE  A ( I » Cl ) OX  OR  G(I) 

= 1 

710 

LET  1=1+1 

720 

IF  I>B  THEN  PRINT  ’HELP1 

730 

LOOP 

740 

LET  Y=A ( I v C ) 

750 

LET  G ( I ) =1 

760 
' 770 

DEFEND 

780 

DEFINE  EXCHANGE 

790 

LET  A ( I f C ) =X 

800 

LET  A ( I r Cl > =Y 

810 

LET  N(X»C)=N(X»C>+1 

820 

LET  N(X»C1)SN(X»C1)-1 

830 

LET  N(Yr  C>=N(YfC)--l 

* 

840 

LET  N<Y»C1)=N(Y»C1)+1 

850 

DEFEND 

855 

860 

DEFINE  CHECK 

862 

rem:  returns  m=i  if  matrix 

NOT  RIGHT 

865 

LET  M=0 

870 

FOR  C=1  TO  K 

880 

FOR  R=1  TO  B 

890 

IF  N(R»CX>1  THEN  LET  M = 1 

900 

NEXT  R 

910 

NEXT  C 

920 

DEFEND 

925 

930 

DEFINE  UNSCRAMBLE 

932 

FOR  J=1  TO  K 

934 

FOR  1=1  TO  B 

936 

LET  A < I » J) =MOD (A(IrJ)-l 

i U)  +1 

938 

NEXT  I 

940 

NEXT  J 

942 

DEFEND 

1000  END 

READY 
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Generalized  Youden  Designs  are  generalizations  of  the  class 
of  two-way  balanced  block  designs  which  include  Latin  squares 
and  Youden  squares.  They  are  used  for  the  same  purposes  and  in 
the  same  way  that  these  classical  designs  are  used,  and  satisfy 
most  of  the  most  common  criteria  of  design  optimality. 

We  explicitly  display  or  give  detailed  instructions  for 
constructing  all  these  designs  within  a practical  range:  when 
v,  the  number  of  treatments,  is/£,25;  and  b-L  and  bg,  the 
dimensions  of  the  design  array,  are  each  £ 50. 
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